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ABSTRACT 

In addition to its great potential for characterizing extra-solar planetary systems the Kepler mis- 
sion is providing unique data on stellar oscillations. A key aspect of Kepler asteroseismology is the 
application to solar-like oscillations of main-sequence stars. As an example we here consider an ini- 
tial analysis of data for three stars in the Kepler field for which planetary transits were known from 
ground-based observations. For one of these, HAT-P-7, we obtain a detailed frequency spectrum and 
hence strong constraints on the stellar properties. The remaining two stars show definite evidence for 
solar-like oscillations, yielding a preliminary estimate of their mean densities. 
Subject headings: stars: fundamental parameters — stars: oscillations — planetary systems 



1. INTRODUCTION 

The main goal of the Kepler mission is to character- 
ize extra-solar planetary systems, p articularly Earth-lik e 
planets in the habitable zone (e.g., iBorucki et all [20091 ). 
The mission detects the presence of planets through the 
minute reduction of the light from a star as a planet 
crosses the line of sight. Several observations of such 
reductions at fixed time intervals for a given star, and 
extensive follow-up observations, are used to verify that 
the effect results from planet transits and to characterize 
the planet. To ensure a reasonable chance of detection 
Kepler observes more than 100,000 stars simultaneously, 
in a fixed field in the Cygnus-Lyra region. Most stars 
are observed at a cadence of 29.4 min, but a subset of 
up to 512 stars can be observed at a short cadence (SC) 
of 58.85 s. Kepler was launched on 6 March 2009 and 
data from the commissioning period and the first month 
of regular observations are now available. 

The very high photometric accurac y required to detect 
planet transits (jBorucki et al.1 120101 : iKoch et all 120101 ) 
also makes the Kepler observations of great interest for 
asteroseismic studies of stellar interiors. In particular, 
the SC data allow investigations of solar-like oscillations 
in main-sequence stars. Apart from the great astrophys- 
ical interest of such investigations they also provide pow- 
erful tools to characterize stars that host planetary sys- 
tems (jKieldsen et al.ll2009D . 

In stars with effective temperature T e g < 7000 K we 
expect to see oscillations similar to tho se observed in the 
Sun (e.g.. lChristcn sen-Dalsgaardll2002| ). excited stochas- 
tically by the near-surface convection. These are acoustic 
modes of high radial order; in main-sequence stars such 
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modes approximately satisfy the asymptotic relation 
v nl ~ Av Q (n + 1/2 + e) - 1(1 + l)D 



(1) 



(|Vandakurovlll967l; lTassou]||1980D . Here u nl is the cyclic 
frequency, n is the radial order of the mode and I is 
the degree, I = corresponding to radial (i.e., spher- 
ically symmetric) oscillations. Also, A^o is essentially 
the inverse sound travel time across the stellar diameter; 
this is closely related to the mean stellar density (p»): 
A^o (p*) 1 / 2 . -D depends sensitively on conditions 
near the center of the star; for stars during the central 
hydrogen burning phase this provides a measure of stellar 
age. Finally, e is determined by conditions near the stel- 
lar surface. This regular form of the frequency spectrum 
simplifies the analysis of the observations, and the close 
relation between the stellar properties and the param- 
eters characterizing the frequencies make them efficient 
diagnostics of the properties of the star. This has been 
demonstrated in the last few years through observations 
of solar-like oscill ations from the ground and from space 
(for reviews, see [Bedding fc Kields"enl 120081 : lAerts et all 
l2lMlG!niiand et al.H20i0alh ~ 

Even observations allowing a determination of Ai/rj 
provide useful constraints on (/?*). With a reliable de- 
termination of individual frequencies (p,) is tightly con- 
strained and an estimate of the stellar age can be ob- 
tained. This can greatly aid the i nterpretation of obser- 
vations of planetary t ransits fe.g.. iGilliland et al"1l2010bl ; 
iNutzman et al.ll2~010f ). We note that photometric obser- 
vations such as those carried out by Kepler are predom- 
inantly sensitive to modes of degree / = 0—2. As indi- 
cated by Eq. ([T|) these are sufficient to obtain information 
about the core properties of the star. 

Ground-based transit observations have identified 
three planetary systems in the Kepler fiel d: TrES-2 
(lO'Donovan et alj 120061: iSozzetti et all 120071), HAT-P-7 
iPal et al.ll2008D. and HAT-P-11 (jDittmann efaLl 12001 
iBakos et al.l 120101 ) . These systems have been observed 
by Kepler in SC mode. Their properties (cf. Table [T]) 
indicate that they should display solar-like oscillations 
at observable amplitudes, and hence they are obvious 
targets for Kepler asteroseismology Here we report the 
results of a preliminary asteroseismic characterization of 
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TABLE 1 
Properties of transiting systems. 



Name 


KIC No 


'1 off (K) 


[Fe/Hj 


L/L 


log(ff) (cgs) 


v sin i 
(kms" 1 ) 


Source 


HAT-P-7 


10666592 


6350 ± 80 


0.26 ± 0.08 


4.9 ± 1.1 


4.07 ±0.06 


3.8 ±0.5 


(a) 






6525 ± 61 


0.31 ±0.07 




4.09 ±0.08 




(b) 


HAT-P-11 


10748390 


4780 ± 50 


0.31 ±0.05 


0.26 ±0.02 


4.59 ±0.03 


1.5 ± 1.5 


(c) 


TrES-2 


11446443 


5850 ± 50 


-0.15 ±0.10 


1.17 ±0.10 


4.4 ±0.1 


2±1 


(d) 






5795 ± 73 


0.06 ± 0.08 




4.30 ±0.13 




(b) 



Note. — Sources: (a): IPal et al.l (2003); (b): lAmmler -von Eif et al.l 1)20091 ) ; (c) : IBakos et al.l l|2010D ; (d) : ISozzetti et al.l (|2007l ) . In some 
cases asymmetric error bars have been symmetrized. 



the central stars in the systems, based on the early Ke- 
pler data. 

2. OBSERVATIONS AND DATA ANALYSIS 

We have analyzed data from Kepler for three 
planet-hosting stars using a pipeline developed for 
fast and robust analysis of all Kepler p-mode dat a 
(|Christensen-Dalsgaard et all 120081 : iHuber et all l2009f ). 
Each time series contains 63324 data points. SC data 
characteristi cs and minor post-pip eline processing are 
discussed in lGilliland et all (|2010cf h In addition a limb- 
darkened transit light curve model fit has been removed 
and 5-cr clipping applied to remove outlying data points 
from each of the time series. The frequency analysis con- 
tains four main steps: 

1. We calculate an oversampled (factor of four) ver- 
sion of the power spectrum by using a least-squares 
fitting. We smoothed the spectrum to 3 /iHz reso- 
lution to remove the fine structure caused by the 
finite mode lifetime. 

2. We correlated the smoothed power spectrum with 
an equally spaced comb of delta functions, sepa- 
rated by Az^o/2, and confined to a Gaussian-shaped 
band with a full width at half maximum of 5Ai/rj. 
We adopted the maximum of this convolution over 
lags between and 0.5 Av as the filter output for 
each Ai/q. 

3. After identifying the peak correlation for the best 
matched model filter and extracting the large sep- 
aration corresponding to this peak we calculate the 
folded spectrum (see Fig. [TJd), i.e., the sum of the 
power as a function of frequency modulo the opti- 
mum large separation (the one corresponding to the 
peak correlation). The summed power is used to 
locate the p-mode structure and identify the ridges 
corresponding to the different mode degrees (based 
on the asymptotic relation). 

4. From the asymptotic relation and the identification 
of mode degrees we finally identify the position of 
the individual p-mode frequencies in the smoothed 
version of the power spectrum; when more than 
one mode is seen near the expected frequency we 
use the power-weighted average of the two peaks. 
Those extracted frequencies and the mode identifi- 
cations are used in the modeling. 
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Fig. 1. — (a) Power spectrum of HAT-P-7 for frequencies between 
300 and 3000 /tHz. The spectrum is smoothed with a gaussian filter 
with a FWHM of 3 /iHz. The noise level at high frequencies corre- 
sponds to 1.1 ppm in amplitude. The white curve is a smoothed 
power spectrum with a gaussian filter (150 /iHz FWHM). A fit to 
the background (dashed white curve) is also shown. The excess 
power and the individual p-modes are evident, (b) Folded power 
spectrum, between 750 and 1500 /tHz, for HAT-P-7 for a large sep- 
aration of 59.22 /tHz. Indicated are the positions corresponding to 
radial modes (I = 0) and non-radial modes with 1 = 1 and 2. The 
measured positions are used to identify the individual oscillation 
modes in panel (a), (c) Bchelle diagram (see text) for frequencies 
of degree 1 = 0, 1, and 2 in HAT-P-7; a frequency separation of 
59.36 /iHz and a starting frequency of 10.8 /iHz were used. The 
filled symbols, coded for degree as indicated, show the observed 
frequencies, while the open symbols are for Model 3 in Table [2] 
minimizing xt- 



For observations with low signal-to-noise ratio it may 
not be possible to identify the individual frequencies. In 
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such cases the analysis is carried through step 2, to de- 
termine the maximum response and hence an estimate of 
the large separation. 
Results on the three individual cases are presented in 

m 

3. MODEL FITTING 

Stellar evolution models and adiabatic oscillation 
frequencies were comp u ted us ing the Aarhus codes 
(|Christensen-Dalsgaardl l2008alTbh . wit h the OPAL 
equation of state dRoeers et al.l [1996) and opacity 
(jlglesias & Rogerslll996fl and the N ACRE nuclear reac- 
tion parameters (|Angulo et al.lll999T) . In some cases (see 
below) diffusion and settling of hel ium were included, 
using the simplified formulation of Mic haud &: P roffitt 
1993). Convection was treated with the iBohm-Vitensa 



1958( 1 mixing-length formulation, with a mixing length 
oml = 2.00 in units of the pressure scale height roughly 
corresponding to a solar calibration. In some models 
with convective cores, overshoot was included over a dis- 
tance of a ov pressure scale heights. Evolution started 
from chemically homogeneous zero-age models. The ini- 
tial abundances by mass Xq and Zq of hydrogen and 
heavy elements were characterized by the assumed value 
of [Fe/H], using as reference a p resent solar surface com- 
position with Z s /X s = 0.0245 JG revessc & Noels 1993) 
and assuming, from galactic chemical evolution, that 
X Q = 0.7679 - 3Z . 

From the observed Av>q, effective temperature and 
composition an initial estimate of the stellar parame- 
ters was obtained using the grid-based SEEK pipeline 
(Quirion et al., in preparation). Smaller grids were then 
computed in the vicinity of these initial parameters, to 
obtain tighter constraints on stellar properties. For HAT- 
P-7 the analysis of the observations yielded frequencies 
of individually identified modes; here the analysis was 
based on 



1 



(obs) 
' i 



(mod) v 

' i 
nl 



N 



1 ^ 

nl 



(2) 



where and are the observed and model fre- 

quencies, a v is the standard error in the observed fre- 
quencies (assumed to be constant) and N is the num- 
ber of observed frequencies. In addition, we considered 
X 2 = xt + Xt' wnere Xt ^ s t ne corresponding normalized 
square difference between the observed and model effec- 
tive temperature. When xt wa s available we minimized 
it along each evolution track and considered the result- 
ing minimum values, and the corresponding value of % 2 , 
as a fun ction of the parameter s characterizing the mod- 
els (see IGilliland et al"1 l2010b. for details). When only 
the large separation A^o could be determined from the 
observations, we identified the model along each track 
which matched Auq and considered the resulting Xr as 
a function of the model parameters. 

4. RESULTS 

4.1. HAT-P-7 

The observed power spectrum for HAT-P-7 is shown 
in Fig.rjju The presence of solar-like p-mode peaks, with 
a maximum power around l.lmHz, is evident. At high 
frequency the noise level in the amplitude spectrum is 



1.1 parts per million (ppm), with some increase at lower 
frequency, likely due to the effects of stellar granulation. 

Carrying out the correlation analysis described in fj2] 
we determined the large separation as A^o = 59.22 /iHz. 
Figure |T|b shows the resulting folded spectrum. This 
clearly shows two closely spaced peaks, identified as cor- 
responding to modes of degree / = and 2, and single 
peak separated from these two by approximately Avq/2, 
corresponding to I = 1. On this basis we finally deter- 
mined the individual frequencies, identifying the modes 
from the asymptotic relation; the final set includes 33 
p-mode frequencies, determined with a standard error 
a v = lAfiUz. These frequencies, corresponding to ra- 
dial orders between 11 and 24, are illustrated in Fig. [TJ: 
in an echelle diagram (see below). 

A grid of models was computed for masses between 
1.41 and 1.61 M Q , [Fe/H] between 0.17 and 0.38, and 
aw = 0,0.1 and 0.2, extending well beyond the end of 
central hydrogen burning. The modeling did not include 
diffusion and settling. At the mass of this star the outer 
convection zone is quite thin, and as a result the set- 
tling timescale is much shorter than the age of the star. 
Including settling, without compensating effects such as 
partial mixing in the radiative region or mass loss, leads 
to a rapid change in the surface composition which is 
inconsistent with the observed [Fe/H]; for simplicity we 
therefore neglected these effects for HAT-P-7. 8 

The computed frequencies wer e corre cted according to 
the procedure of iKieldsen et alj (|2008f ) for errors in the 
modeling of the near-surface layers, by adding a(v/v^) h 
where a = 0.1158 /iHz, v = 1000 /iHz and b = 4.9. As 
discussed in Sj3[ for each evolution track, characterized 
by a set of model parameters, we minimized the depar- 
ture xt °f the model frequencies from the observations, 
defining the best model for this set. 

We first consider xt as a function of the effective 
temperature of the models (Fig. EH). It is evident 
that there is a clear minimum in xt\ this is consi stent 
with the determination of T e g by IPal et all (|2008l ) but 
not with the somewha t highe r temperature obtained by 
lAmmler-von Eif et alJ (|2009h (see also Table [TJ. Thus 
in the fo l lowing we use the observed quantities from 
IPal et alJ pOOaf . 

Since the frequencies to leading order are determined 
by the mean stellar density (/?*), Fig. [2p,c show xt an d 
X 2 as functions of (p*). It is evident that the best-fitting 
models occupy a narrow range of (p*), with a well-defined 
minimum. Fitting a parabola to x 2 m panel (c) we obtain 
the estimate (/?»} = 0.2712 ± 0.0032gc m - 1 . In Fig. EH 
X 2 is shown against model age. Here the variation with 
model parameters is substantially stronger, resulting in 
a greater spread in the inferred age; in particular, it is 
evident, not surprisingly, that the results depend on the 
extent of convective overshoot. From the figure we esti- 
mate that the age of HAT-P-7 is between 1.4 and 2.3 Gyr. 

Examples of evolution tracks are shown in Fig. (3[ pa- 
rameters for these models are provided in Table EJ They 
were chosen to give the smallest xt f° r each of the three 
values of a ov considered. Also shown are the locations 

8 Artificially suppressing settling in the outer layers, while in- 
cluding diffusion and settling in the core, leads to results that are 
very similar to those presented here. 
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TABLE 2 

Stellar evolution models fitting the observed frequencies for HAT-P-7. 



No 


M„/M 


Age 


Zo 


X 


a ov 


R,/Rq 


<P*> 


Teff 


WL© 


xi 


x 1 






(Gyr) 










(gem 3 ) 


(K) 








1 


1.53 


1.758 


0.0270 


0.6870 


0.0 


1.994 


0.2718 


6379 


5.91 


1.08 


1.21 


2 


1.52 


1.875 


0.0290 


0.6809 


0.1 


1.992 


0.2708 


6355 


5.81 


1.04 


1.04 


3 


1.50 


2.009 


0.0270 


0.6870 


0.2 


1.981 


0.2718 


6389 


5.87 


1.00 


1.24 



Note. — Models minimizing xi (cf- Eq.[2]| along the evolution tracks, illustrated in Fig. [3] The models have been selected as providing 
the smallest xi f° r each of the three values of the overshoot parameter a ov . The smallest value of xi ls obtained for Model 3. 
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Fig. 2. — Results of fitting the observed frequencies to a grid of 
stellar models (see text for details). Plusses, stars and diamonds 
correspond to models with a ov = (no overshoot), 0.1, and 0.2. (a) 
Minimum mean square deviation xi °f the- frequencies (cf. Eq. [2j 
along each evolution track, against the effective temperature T ff 
of the corresponding models. The vertical da shed and dotted lines 
indicate the effective temper atures found by Pal et al. ( 20081) and 
Ammlcr-von Eif et al. (2009). (b) Minimum mean square devia- 
tion xi against the mean density (p„) of the corresponding models, 
(c) As (b), but showing the combined x' 2 ■ (d) X 2 against the age 
for the models that minimize xi I the different ridges correspond to 
the different masses in the grid, the more massive models resulting 
in a lower estimate of the age. 




7000 6800 6600 6400 6200 6000 

T,„ (K) 

Fig. 3. — Theoretical HR diagram with selected evolutionary 
tracks, corresponding to the models defined in Table 2. The in- 
dicate the models along the full set of evolutionary sequences mini- 
mizing the difference between the computed and observe d frequen- 
cies. T he box is centered on the L and T ef f as given bv lPal et alj 
(2008), with a size matching the errors on these quantities. 

of the models minimizing xt along each of the computed 
tracks; these evidently fall close to a line in the HR di- 
agram, corresponding to t he small range i n (p*). The 
range of luminosities, from I Pal et al.l (|2008f ). is based on 
modeling and hence has not been used in our fit; even 
so, it is gratifying that the present models are essentially 
consistent with these values. Also, as indicated by Fig. [5^ 
and Table[2j the be st-fitting models are close to the value 
of T eff obtained bv lPal et all (|2008h . 

The match of the best-fitting model (Model 3 of Ta- 
ble [2|) to the observed frequencies is illu strated in a so- 
called echelle diagram (Gr ec et al.| [T983) in Fig. [TJ:. In 
accordance with Eq. ([T| the frequency spectrum is di- 
vided into slices of length Av, starting at a frequency of 
10.8 //Hz; the figure shows the location of the observed 
(filled symbols) and computed (open symbols) frequen- 
cies within each slice, against the starting frequency of 
the slice; the model results extend to the acoustical cut- 
off frequency, 1930 pHz, of the model. There is clearly a 
very good overall agreement between model and obser- 
vations, including the detailed variation with frequency 
which reflects the frequency dependence of the large sep- 
aration, as a possible diagnostics of the outer layers of 
the star (e.g., iHoudek fc Gough|[2007l) . 

We have finally made a fi t of the inferred (p*), as 
well as T cB and [Fe/H] from iPal et"al] (|2008l ). to com- 
puted evolutionary tr acks from the Yonsei-Yale compi- 
lation (|Yi et alJ 120011 ) . This was based on a Markov 
Chain Monte Carlo analysis to obtain the statistical 
properties of the inferred quantities (see iBrown I I2010L 
for details). This resulted in M = 1.520 ± 0.036 M , 
R = 1.991 ± 0.018 R Q and an age of 2.14 ± 0.26 Gyr. 
We note that the age estimate reflects the specific as- 
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sumptions in the Yonsei-Yale evolution calculations; as 
indicated by Fig. [2ji the true uncertainty in the age de- 
termination is likely somewhat larger. 

4.2. HAT-P-11 

For HAT-P-11 the oscillation amplitudes were much 
smaller than in HAT-P-7, as expected from the general 
scalin g of amplitudes with stel lar mass and luminosity 
(e.g. JKieldsen & Bedchnell 19951 ). Thus with the present 
short run of data it has only been possible to determine 
the large separation A^o = 180.1 /iHz from the maximum 
in the correlation analysis. We have matched this to a 
grid of models, including diffusion and settling of helium, 
with masses between 0.7 and 0.9 M© and [Fc/H] between 
0.21 and 0.41. These models provide a good fit to the 
observed T e g and L/L©; note that in the present case the 
luminosity is based on a reasonably well-determined par- 
allax. We have determined an estimate of (p» } by averag- 
ing the results of those models which match the observed 
A^o and lie within 2 standar d deviations (±100 K) from 
the value of T e g provided bv lBakos et al.l (|2010f ); the re- 
sult is (p*) = 2.5127 ±0.0009gcm" 3 . Although the for- 
mal error is extremely small, owing to a tight relation 
between the large separation and the mean density for 
stars in this region in the HR diagram, the true error is 
undoubtedly substantially larger. In particular, we ne- 
glected the error in the determination of A^o and these 
data have not allowed a correction for the systematic 
errors in the modeling of the near-surface layers of the 
star. 

4.3. TrES-2 

Here also we were unable to determine individual fre- 
quencies from the present set of data. The expected am- 
plitudes are smaller than for HAT-P-7, and the noise 
level higher due to the fainter magnitude of TrES-2. 
The correlation analysis yielded two possible values of 
Av n : 97.7/iHz and 130.7/uHz. For this star (p,) has 
been d etermined from the a nalysis of the transit light 
curve. ISozzetti et all (|2007l) obtained (p*) = 1.375 ± 
0.065gcm" 3 , while iSouthworthl (|2009h found <p«) = 
1.42 ± 0.13 gem" 3 . From the scaling with (p*) 1 / 2 the 
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smaller of the two possible values of A^o is clearly incon- 
sistent with these values of (p*), while Az/o = 130.7 /iHz 
yields models that are consistent w ith the observed T e g 
and log(<7) of ISozzetti et al.l (|2007l ) as well as with these 
values of the mean density. Here we considered a grid 
of models with helium diffusion and settling, masses be- 
tween 0.85 and 1.1 M© and [Fe/H] between -0.25 and 
—0.05. Determining again the mean value of (p*) for 
those models that matched Ais an d had T e g within two 
standard deviations of the value of ISozzetti et al.l (|2007f ) 
we obtained (p») = 1.3233 ± 0.0027gcm" 3 . As in the 
case of HAT-P-11 the true error is likely substantially 
higher. 

5. DISCUSSION AND CONCLUSION 

The present preliminary analysis provides a striking 
demonstration of the potential of Kepler asteroseismol- 
ogy and its supporting role in the analysis of planet hosts. 
These stars will undoubtedly be observed throughout the 
mission and hence the quality of the data will increase 
substantially. For HAT-P-7 the detected frequencies are 
already close to what will be required for a detailed anal- 
ysis of the stellar interior, beyond the determination of 
the basic parameters of the star. Thus here we can look 
forward to a test of the assumptions of the stellar mod- 
eling; the resulting improvements will further constrain 
the overall properties of the star, in particular its age. 
Also, given the observed v sini we expect a rotational 
splitting comparable to that observed in the Sun, and 
hence likely detectable with a few months of observa- 
tions. For the other two stars there is strong evidence 
for the presence of solar-like oscillations; thus continued 
observations will very likely result in the determination 
of individual frequencies and hence further constraints 
on the properties of the stars. 
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NASA's Science Mission Directorate. We are very grate- 
ful to the entire Kepler team, whose efforts have led to 
this exceptional mission. The present work was sup- 
ported by the Danish Natural Science Research Council. 
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